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D 
^n : Abstract 

O ' Analytical calculations of the dynamics of a curved domain wall in a nematic 



liquid crystal are performed. The core of the wall is assumed to form a 

cylinder, whose axis coincides with the direction of an external magnetic field. 

^' 

Ch ' The equation of motion for the nematic director field is solved in a comovine; 

O" 

coordinate frame by applying a polynomial expansion of the tilt angle with 
K^ respect to the radial distance from the wall core. Starting from a cylindrical 

ly-^ ' domain wall at rest as initial conditions, the shrinking of the cylinder and 

o' 

CN ' the change of the wall width is analysed in detail. In particular, we find that 



the Neel wall decays faster than the Bloch wall, in agreement with energy 
C^ ' considerations. 



o 
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I. INTRODUCTION 

It is a well-known fact that nematic liquid crystals reveal a large variety of defect struc- 
tures. Such defects are very interesting from the theoretical point of view ^^ . Moreover, 
they are important for applications, e.g., they play an crucial role in the switching mecha- 
nisms in modern liquid crystal display devices (surface-stabilized ferroelectric liquid crystal 
cells) IQ-0 . Therefore, a quantitative investigation of defect dynamics is desirable. 

The types of defects in the nematic director field cover zero- and one-dimensional objects 
(point and line defects, respectively) that arise due to the presence of impurities or surfaces 
IPJ^. Two-dimensional sheets, i.e. walls, are unstable in an isolated nematic. They can 
be stabilized, however, by imposing an external magnetic field |^. The situation is similar 
to domain walls in a ferromagnetic material where regions of different spin orientations 
are separated by a thin boundary layer |T0|. In a nematic, the director can align parallel 



or antiparallel to the external magnetic field. Due to fact that the director "lacks the 
arrowhead", these two situations are energetically equivalent. It may happen, that two 
spatial regions of antiparallel director orientation will form, with a domain wall between 
them p. 

Domain walls in nematics due to the existence of an external field are static in a planar 
geometry [^ . A qualitatively new behaviour occurs, when the shape of the wall is curved. In 
this case the wall becomes a dynamical object, i.e. its core starts to move, whereas the wall 
width is changing. The present paper is devoted to a quantitative study of the dynamics 
of domain walls in a cylindrical geometry. We consider the case of director reorientation 
through a Bloch and Neel wall ||Tl],|l2|. From a methodological point of view, in our calcu- 



lations we use the so-called polynomial approximation - it has turned out to be very fruitful 



in studies of defect dynamics in relativistic scalar field theories [jT3[. Due to the cylindrical 
symmetry, we can find explicit analytical solutions to equations of motion obtained within 
that approximation. 

The organization of the paper is as follows. In section II the equation of motion for the 



director field is derived on the basis of tlie Ericksen-Leslie tlieory p^ , [T5[] for the case of a 
cyhndrical domain wall of Bloch or Neel type in a nematic liquid crystal. Section III develops 
the polynomial expansion as an analytical method to solve approximately the dynamics of 
the director field. Within this approach, we obtain closed expressions for the dynamics of 
the wall core and width. In section IV some selected results together with their discussion 
are presented. Section V contains an outlook on possible extensions of our investigations. 

II. DIRECTOR EQUATIONS OF MOTION FOR CYLINDRICAL DOMAIN 

WALLS 

A. Comoving polar coordinates 

For all calculations we restrict ourselves to the cylindrical geometry as pointed out in 
Figs. 1 and 2. The cylinder axis is determined by the orientation of the external magnetic 
field H which coincides with the z-axis of our coordinate frame, i.e. H = Hq z. {z denotes 
the unit vector along z.) As already mentioned above, a domain wall occurs for parallel and 
antiparallel director orientation with respect to the external field in different spatial regions. 
In the case of axial geometry, we assume the director to be aligned along the positive or 
negative 2;- axis {up, down) inside or outside the cylinder, respectively. Now the director 
can reorient from up to down in either (p — 2;— plane or p — z— plane, where (p, 0, z) are 
cylindrical coordinates. We refer to these two cases as Bloch (Fig. 1) and Neel (Fig. 2) wall, 
taking over the nomenclature from planar domain walls in ferromagnetics [1^ and nematics 



For our calculations it will be advantageous to work within a comoving cylindrical coor- 
dinate frame. The radial coordinate consists of the wall radius p{t) and the radial distance 
^ from the wall. These are related to cartesian coordinates according to 

^ = (pit) + COS 0, y= (pit) + sin 0. (1) 

The Nabla operator, expressed in the comoving coordinates, reads 



The "spatial" time derivative d /dt is related to the "material" {i.e. co moving) time deriva- 
tive D /Dt by the chain rule between the coordinate sets Xi = {x, y, z) and ^q, = (,^, 0, z) 
(laboratory and comoving coordinates): 



d D dxi 
di^Dt~ ~dt 



Co 



dxi d^, 
The relation above is simplified for cylindrical geometry 



9^ d ^3^ 



- = — -•— (A) 

dt Dt ^ di' ^ ^ 

For the Bloch and Neel wall we can describe the director orientation by the tilt angle field 
Q{C,,t) which measures the angle between the local director and the z-axis. In terms of the 
comoving frame, the director is in the — 1;— plane for the Bloch wall and in the ^ — z— plane 
for the Neel wall: 

riB = sin e(^, t) ^(0, t) + cos e(^, t) z, (Bloch wall) (5) 

njv = sine(^,t)^(0,t) + cose(^,t)z. (Neel wall) (6) 

B. Free energy density 

In the framework of the Ericksen-Leslie theory [0,^ the dynamics of the director field 
follows from the balance between elastic, magnetic and viscous torques. The latter are 
determined by the temporal change in the director, whereas elastic and magnetic torques 
are obtained as a variational derivative from a free energy density J-'. Thus the equation of 
motion for the director components generally reads 

dui 5J^ _ dT dT . . 

dt 6ni '' d^djUi) drii ' 

In Eq. (|^) 7i is the rotational viscosity of the liquid crystal. In fact, Eq. (0) should contain 
an additional term An^, with the Lagrange multiplier A, accounting for the constraint of 



n being a unit vector. In the cylindrical geometry under consideration, however, it turns 
out, that this constraint can be incorporated much more easily by starting from (J^) {i.e. 
without the Lagrange multiplier). Later on we shall eliminate the variation of the length of 
the director by taking appropriate projections, as will be demonstrated below. 

The free energy density JF consists of an elastic and a magnetic part (JF = J^eiast + ^mag) ■ 
For the elastic free energy density we take the Oseen-Zocher- Frank expression ||l6Hl8|, that 
contains splay, twist and bend deformations of the director field. 

J^eiast = - Kn (divn)^ + - K22 (n ■ curl n)^ + - K33 (n x curl nf. (8) 

22 2 

In (|^) Kii, K22 and -^"33 denote elastic constants of the nematic liquid crystal. 

The magnetic free energy density couples the director n to the magnetic field H via the 
anisotropy of the magnetic susceptibility Ax (a^o means the magnetic field constant). 

^mag = — /ioAx(n-if)^ (9) 

For a cylindrical Bloch or Neel wall we can express the free energy density explicitly in 
comoving coordinates, 

^ 1 ^^ /^, cos 6 sin6\ 1 ^^ sin^ 9 

^B = -i^22 e'+ ——] +-K33 



\ P+e ; 2 ''{p+0' 

-- /io Ax H^ cos^ e, (Bloch wall) (10) 



2 

2 

'2 „• 2 , 



1 - - / „ , „ sin B \ 1 



jr^ = - Ku e' cose + + - K33 sin^ 9 

2 V p + ^J 2 

-- /io Ax H^ cos^ e. (Neel wall) (11) 



In (|T0|) and (|Tl|) O' stands for the radial derivative dQ/d^. 



C. Director equation of motion 

Returning to Eq. (|^ we now have to account for the constraint of the director normal- 
ization (|n| = 1). To this aim we write the left-hand side (LHS) of (0) explicitly in comoving 
coordinates. 



LHSb = -71 (6 - p e') (^ X n), (Bloch wall) (12) 

LHSjv = 7i (e - p e') (0 X n). (Neel wall) (13) 

Here denotes the material time derivative DQ/Dt. 

For further considerations we express the right-hand side (RHS) of (|^) in cartesian co- 
ordinates. Its components read 

RHSj = Ku di (divn) — K22 (n ■ curln) (curln)j 
—K22 {curl [(n ■ curln) n]}j 
-K33 {{uk dkUj) diUj - dk[nk{nj djUi)]} 
+fioAx{n-H)Hi. (14) 

The constraint of the director normalization |n| = 1 means that we have to discard the 
(infinitesimal) variation of the length of the director 6n = e(r) n. (e(r) is a small number 
everywhere.) The corresponding variation of the total free energy F (which is the volume 
integral over the free energy density JF) reads 



SF = / dV -— ■ 5n = / d'r e(r) n ■ -— = 0. (15) 

J on J on 

From Eq. ( |15|) it is obvious that the variation of the director length is related to the 
projection of the variational derivative of the free energy density [i.e. of the RHS ([T^)) 
onto the director. Moreover, from (0), (|l^) we recognize that the LHS, projected onto the 
director, yields zero. Thus we conclude that by discarding the projection of Eq. (|^) onto the 
director we properly take into account the director constraint! 

Additionally, after some lengthy calculations we find that for the Bloch and Neel wall 
the projection onto ^ and 0, respectively, yield identically zero. Thus the scalar equation 
of motion for the tilt angle field 0(.C, t) can be obtained by projecting (|^) onto the the third 
linearly independent direction, which is {$, x n) for the Bloch wall and (0 x n) for the 
Neel wall. After performing the projections and changing for the comoving coordinates, we 
obtain the equations of motion for the tilt angle field. 
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7i(e-pe') = K22e" + ir22 ® 



sin 6 cos 6 

--f^22 



-2 (7^33 - K: 



iP + 0' 

cos O sin^ O 



22 J 



(p + 



2 



-/io Ax ^0 sin e COS e, (Bloch wall) (16) 

7i(e-peo = Kncose^|-^^[(p + Osine]| 

sin 9 cos 9 

-/io Ax i^o sin 9 cos 9. (Neel wall) (17) 

In agreement with the pictorial visualization of Figs. 1 and 2 the cylindrical Bloch wall 
contains no splay deformations, whereas in the Neel wall twist deformations are absent. 



III. SOLUTION FOR THE DIRECTOR DYNAMICS 
A. Polynomial expansion 

The method to construct an approximate solution of the equations of motion (|TBp, (|T^ 
is the polynomial expansion of the field which has been developed in previous papers ||13[| . 
We shall specify it here to the case of cylindrical domain walls in nematic liquid crystals. 

The key idea is to take the spatial dependence of the tilt angle field 9(^,t) as a Taylor- 
like expansion with respect to the wall distance ^ around its "core" value 9(^ = 0, t). The 
wall core corresponds to the radius of the cylinder {i.e. ^ = 0), where the director is oriented 
perpendicular to the 2;-axis, thus 9(^ = 0, t) = | Vt. The temporal evolution of the field is 
governed by the expansion coefficients. 

9(e, t) = - + a{t) e + - b{t) e + - c{t) e- (18) 

2 2 6 

The polynomial expansion cannot be truncated at arbitrary order, because we have to glue 
the bulk solution for 9(.^,t) smoothly to the boundary conditions. This is not always 



possible. However, it has been demonstrated for the case of cyhndrical domain walls in 
a relativistic field theory (where the same symmetry considerations are valid as in our 
case), that the third-order polynomial expansion (p!8|) guarantees the compatibility of the 



bulk solution with the boundary conditions [13|. Of course, for higher order of the Taylor 



expansion the accuracy of the calculations will increase. 



From ([T8|) we can easily obtain the spatial and temporal derivatives of the tilt angle field 
Q{C,,t). The trigonometric and fractional expressions that occur in the equations of motion 
(|16|), (\[7\} are expanded as well up to third order in the wall distance ^. All these expansions 
are then inserted into (|16D and ([171) . By comparison of the coefficients of subsequent orders 
in ^ we obtain a set of ordinary differential equations for the expansion coefficients a{t), b{t) 
and c{t). For the third-order polynomial expansion ([18]) such a comparison yields meaningful 
results for zeroth and first order in ^. Thus from our calculations we obtain two equations 



that read as follows: 



7i p a = -K22 b-K22-, (Bloch wall) (19) 

71 (d - p6) = K22 C + K22- + 2 (7^33 - K22) 4 

P P 

+Po Ax H^ a, (Bloch wall) (20) 

71 p a = -7^33 6 - 7^33 -, (Neel wall) (21) 

b 
71 (d - p 6) = 7^33 c + K33 - 

P 

-{K33-Ku)^-{K33~Kn)a' 
P^ 

+Po Ax Hq a. (Neel wall) (22) 



B. Boundary conditions 

Next we have to incorporate the boundary conditions. These are determined by the 
director alignment parallel and antiparallel to the external magnetic field inside and outside 
the wall, respectively. The (instantaneous) wall thickness is given as ^o(^) + 'Ci(i)? where 
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^o(^) measures the outward distance from the core to the outer edge of the wall, whereas 
C,i{t) means the core-to-edge distance towards the cylinder axis. Therefore the boundary 
conditions read 

e{^o,t) = n, 0(-a,t) = O, (23) 

e'(eo,t) = 0, 0'(-^i,t) = O. (24) 

Eqs. ( p3D determine the director orientation at the edge of the wall, and Eqs. (0) are the 
conditions for the smoothness of the solution. 

We can immediately express the boundary conditions in terms of the polynomial expan- 
sion for the tilt angle field (|TBp: 

a^o + lbeo + -^ceo = -^, (25) 

a + b^o + -ceo = 0, (26) 

2 

-a^^ + -be,--ce, = --, (27) 

2 6 2 

a-b^i + -c^l = 0. (28) 

2 

The equations above form a set of inhomogeneous linear equations for the expansion co- 
efficients a(t), b(t), c(t). The solubility conditions yield relations between the expansion 
coefficients and the wall partial widths ^o(^) and ^i(t). 

(29) 
(30) 

We note that the partial widths of the wall have turned out identical, and that there is no 
quadratic term in the polynomial expansion. 



^0 = 6, 


"=4&' 






37r 


32 n 


6 = 0, 






2?g 


97r2 



These relations must hold also in the bulk equations ([T9|) -(22) in order to fulfill the 



boundary conditions. Introducing them, ([T9|)-(p2D become simplified. 

71 p = -:^, (Bloch wall) (31) 

32 n , , a 

71 « = r ^22 0^ + 2 iKs3 - K22) ^ 





+yUo Ax Hi a, (Bloch wall) 




(32) 


7iP = 


^^ , (Neel wall) 




(33) 


7i« = 


= -{— K33 + K33 - Kn) a' - (i^33 - 

97r2 


- ^11) — 
P 






+/io Ax Hi a. (Neel wall) 




(34) 



Equations ( pi] ) and ( P3| ) describe the time evolution of the wall core, i. e. the cylinder radius, 
whereas (|3^) and (^) govern the dynamics of the wall half-width, for the Bloch and Neel 
wall, respectively. 

C. Solution for the w^all dynamics 



The coupled equations of motion (|3TD-(p^) can be solved analytically. First we solve 
(pip and (|33D for the dynamics of the wall core. It obeys a square root law: 

(35) 

(36) 

V 7i 

Qualitatively, starting from a wall radius po, the cylindrical domain wall will shrink. 



p(t) = = 


^P§-2^, 


(Bloch wall) 


Pit) = = 


J Pi i^^^t. 


(Neel wall) 



This behaviour is similar to what has been found for relativistic field theories 131. For usual 



nematics 7^22 is the smallest of the elastic constants. Therefore the (|35D and (|36|) indicate 
that the decay time of the Bloch wall is larger than that of the Neel wall. 

Now we can use the solutions (p5D, (|36| ) to eliminate the p-dependence from the equations 
of motion for the wall width (P^D, (0). The latter then read 

(37) 



a = A{t)a- Ba^ 



where the abbreviations stand for 



A^t)^nK..-j<..)^H^jj, 



B 



71 Po 

32 K' 



2 Root 



71 



0' 



22 



97r2 7i 



(Bloch wall) 



(Bloch wall) 



(38) 
(39) 
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Ait) 
B 



(i^33-irn) , /ioAx 



7i pg - 2 i^33 1 7i 

32 7^33 i^33 - -^11 



+ '-^^-^ Hi, (Neel wall) 



(Neel wall) 



(40) 
(41) 



97r2 7i 7i 

Eq. (|37|) is an ordinary differential equation of Bernoulli type. By the substitution a = l/va 
it can be transformed into an inhomogeneous linear differential equation 



a + 2A{t)a = 2B, 



whose general solution is obtained from a two-fold integration: 



(42) 



~a{t) 



M{t) 



[2BN{t) + ao] 



M{t) =exp(/ dt'2A(t')), 

^0 



(43) 



(44) 



N{t) = / dt'M(t'). (45) 

Jo 

fio is the initial value for the auxiliary variable a{t). This variable is related closely to the 
wall half- width ^o(^) (Eq-(51) below). 

For our case of the cylindrical domain walls all integrals can be solved analytically. 



M{t) =e"*(l-/?t) 
N(t) "- '^ 



a/ 13 / o\ S+1 



P \q;/ 
where further abbreviations are introduced 



1(6 + 1,^) -^{5 + 1,^ -at) 



(46) 
(47) 



a 



'^^^Hl (3=^, 5 = -2^23 +2, (Blochwall) 



7iPo 



(48) 



2/ioAx TTl n_ 2 7^33 

— ^^ — -"o, P — ,, '1- 



7iPo 



5 = 1- 






33 



(Neel wall) 



(49) 



In ( ^ ) 7(1/, x) stands for the incomplete Gamma function, whose integral representation is 



7(1^,0;) = / dte 't 

JO 



-t ,u-l 



(50) 
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From the solution for the auxihary variable d{t) the time evolution of the wall half-width 
C,o{t) is obtained by simple back substitution, 

Stt Sir 



Ut) 



4a(t) 



4 



~ait). 



(51) 



IV. QUANTITATIVE RESULTS 

The analytical results for the dynamics of cylindrical domain walls in nematics that 
have been obtained in the previous section can be specified to real materials. For further 
considerations we choose the parameters that enter the solutions according to the nematic 
phase of PAA (p-azoxyanisole) at 120°C [|,|1 (Table I). 



splay elastic constant Kn 


7.0 • 10-12 N 


twist elastic constant K22 


4.3 ■ 10-12 N 


bend elastic constant ^33 


1.7 -10-" N 


rotational viscosity 71 


6.7-10-3 Nm/s 


magnetic anisotropy /iq Ax 


1.21-10-^ 



TABLE I: Material constants of PAA at 120°C. 



The magnetic field strength Hq has been chosen 500 Oersted, according to a magnetic flux 
density Bq = iiqHq = 0.05 T. 

The initial configuration is a cylindrical domain wall at rest. At zero time the cylinder 
radius is by two or three orders of magnitude larger than the wall half-width: po = p{t = 
0) = 0.1mm, with ^0 = ^o(t = 0) = 1/im or 0.1 /im. 



Fig. 3 shows the shrinking of the wall according to the square- root law (|35|) , (0). The 
actual decay time is the moment when the cylinder radius touches the edge of the wall. 
However, from the temporal evolution of the wall width (see below) it is obvious that the 
decay time is only slightly overestimated when taking the time for which the cylinder has 
shrinked to zero. 
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2 
Tdecay = ^, (Bloch Wall) (52) 

2 
Tdecay = ^. (Ned Wall) (53) 

The main feature of the Bloch wall are twist deformations, whereas the Neel wall rather 
consists of bend deformations. The decay times (p^, (^) depend on the elastic constants 
i^33 and K22'- the Neel wall decays about four times faster than the Bloch wall. 

Interestingly, the dynamics of the wall width reveals two separated processes on different 
time scales. First a rapid change occurs from the initial half-width to a metastable state 
that does not depend on the initial condition (Figs. 4 and 5). The metastable half-width can 
be related to the magnetic coherence length ^mag- For planar geometry the director would 



reorient in an external magnetic field by a twist deformation on a length scale [jT9[ of 



Uag = \ — T— -fT (54) 

Introducing the parameters given at the beginning of this section, we obtain ^^ag = 0.15 /im 
for the half-width of a planar Bloch wall. However, according to Figs. 4 and 5 the metastable 
value for the half-width of the cylindrical Bloch wall is .^Bioch = 0.21 yum. The discrepancy 
is due to the curvature that gives rise to additional bend deformations. Assuming that a 
similar law as Eq. (|54D also holds for the curved wall, we can extract an effective elastic 
constant -ft'eioch = 8.4 ■ 10~^^ N. For the curved Neel wall the coherence length (Figs. 4 and 
5) is ,^Neei = 0.68 /xm, corresponding to an effective elastic constant of i^Neei = 8.85 ■ 10^^^ N. 
This means that the elastic energy content of cylindrical domain walls is about one order of 
magnitude larger for the Neel geometry as compared to the Bloch-like reorientation. 

As already stated above the effective coherence length ^Bioch or ^nccI denotes a metastable 
state. After some seconds, when the cylinder has shrinked so that its radius p(t) is of the 
order ^o(^)) our solution for the wall half-width breaks down. Formally, our expressions for 
the Bloch wall reveal an implosion (Fig. 6), whereas for the Neel wall we obtain an explosion 
(Fig. 7) of the wall width. This is due to the opposite sign of the quantity 5 in ( |I5| ) and 
(13) and thus depending on the relative magnitude of the elastic constants (c/. (^Sj) , (^9D ). 

13 



Actually, one should remember that for wall half-widths larger than the cylinder radius 
the boundary conditions (|2^ ) are not correct, and then our solutions loose their physical 
meaning. Thus, within the framework of the polynomial expansion we are able to study 
the decay of the cylindrical domain wall until its radius becomes approximately equal to its 
half-width. 

V. REMARKS 

1. Our solutions for p{t) and ,^o(^) give a rather detailed description of the time evolution 
of cylindrical domain walls. First of all, it would be interesting to compare our theoretical 
predictions with experimental results. Furthermore, because our expressions explicitly show 
how the dynamics of the domain wall depends on the material constants Kn, K22, -ft'sa, 
fio and 7i, they could be useful for the determination of these constants by observing the 
shrinking of the cylindrical domain walls. 

2. It is possible to generalize our calculations to domain walls of a more general shapes, 
e.g. we could allow for a modulation along the 2;-axis or for non-circular sections by the x- 
y-plane. For the relativistic scalar field theory this has been performed in the second of the 
references |13| . The calculations in the general case the are much more cumbersome than for 
strictly cylindrical symmetry, while the main ideas of the polynomial approximation remain 
unchanged. For this reason we have chosen to restrict our presentation to the cylindrical 
geometry. 
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FIGURES 
FIG. 1. Geometry and coordinates for a cylindrical Bloch wall. Projection onto a plane per- 
pendicular to the cylinder axis. 

FIG. 2. Geometry and coordinates for a cylindrical Neel wall. Projection onto a plane perpen- 
dicular to the cylinder axis. 

FIG. 3. Temporal evolution of the radius of a cylindrical domain wall in PAA at 120°C. Initial 
wall radius p{t = 0) = 0.1mm. Solid: Bloch wall, dashed: Neel wall. 

FIG. 4. Temporal evolution of the half-width of a cylindrical domain wall in PAA at 120°C. 
Initial wall half-width p{t = 0) = 1 pm. Solid: Bloch wall, dashed: Neel wall. 

FIG. 5. Temporal evolution of the half-width of a cylindrical domain wall in PAA at 120°C. 
Initial wall half-width p{t = 0) = 0.1 ;um. Solid: Bloch wall, dashed: Neel wall. 

FIG. 6. Catastrophic behaviour of the width of a cylindrical Bloch wall in PAA at 120°C. 

FIG. 7. Catastrophic behaviour of the width of a cylindrical Neel wall in PAA at 120°C. 
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0.8 - 



0.6 



0.4 - 



0.2 












0.05 



0.1 



0.15 



0.2 



time 



ms 



J. Stelzer, H. Arodz: Figure 5 
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J. Stelzer, H. Arodz: Figure 6 
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J. Stelzer, H. Arodz: Figure 7 
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